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We study mesoscopic fluctuations and weak localization correction to the supercurrent in Joseph- 
son junctions with coherent diffusive electron dynamics in the normal part. Two kinds of junctions 
are considered: a chaotic dot coupled to superconductors by tunnel barriers and a diffusive junction 
with transparent normal-superconducting interfaces. The amplitude of current fluctuations and the 
weak localization correction to the average current are calculated as functions of the ratio between 
the superconducting gap and the electron dwell energy, temperature, and superconducting phase 
difference across the junction. Technically, fluctuations on top of the spatially inhomogeneous prox- 
imity effect in the normal region are described by the replicated version of the cr-model. For the 
case of diffusive junctions with transparent interfaces, the magnitude of mesoscopic fluctuations of 
the critical current appears to be nearly 3 times larger than the prediction of the previous theory 
which did not take the proximity effect into account. 
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I. INTRODUCTION 

At low temperatures, conductance of metals is due 
to electron scattering on impurities. The wave na- 
ture of electron motion reveals in a number of quan- 
tum interference effects: the weak localization (WL) 
correction^ to the classical Drude conductance, and uni- 
versal sample-to-sample conductance fluctuations.^'^ In 
mesoscopic samples whose size does not exceed the phase 
coherence length, L^{T), the magnitude of conductance 
fluctuations characterized by the root mean square (rms) 
6G is independent on the system size and degree of dis- 
order and is of the order of the conductance quantum 
Gq = /irh. Weak localization corrections and univer- 
sal conductance fluctuations have attracted considerable 
interest since the 80s' both from the experimental and 
theoretical sides.'' 

Two years after discovery of conductance fluctuations 
in metals, Altsliuler and Spivak^ applied the same idea 
to fluctuations of the supercurrent in Josephson junc- 
tions formed of a diffusive normal metal (N) placed be- 
tween two superconducting (S) leads. They considered 
the limit of long junctions, when the Thouless energy 
Et = hD/L^ {D is the diffusion constant and L is the 
length of the junction) is much smaller than the super- 
conducting gap A in the leads. In particular, it was found 
in Ref. 5, that for quasi-one-dimensional long wires, 
mesoscopic fluctuations of t he criti cal current are charac- 
terized by the rms 5Ic = ^/3C{3) eET/nh = OMeEr/h 
at zero temperature. 

The theory of Ref. 5 was based on the standard dia- 
grammatic technique operating with soft diffusive modes 
- diffusons and Cooperons - but Andreev reflection at the 
NS interface^ was described by the linear phenomenolog- 
ical boundary conditions on diffusive modes^. However, 
the proximity effect in SNS systems is known to be es- 
sentially nonperturbative at low energies, where Cooper 
pairs penetrating from a superconductor strongly mod- 



ify electronic properties and open a minigap E^ in the 
normal region.'* This strong perturbation of the metal- 
lic state can be described only with the help of the full 
set of nonlinear Usadel equations^ (in diffusive systems) . 
Thus, the treatment of Ref. 5, which effectively consid- 
ered mesoscopic fluctuations of the supercurrent on top 
of a uniform metallic state without a minigap, is not ac- 
curate and should be reconsidered taking the proximity 
effect into account nonperturbatively.^" 

Extensive studies of the proximity effect in SNS 
systems^ demonstrate that the critical Josephson cur- 
rent is related to the minigap induced in the normal re- 
gion: Ic GEt,/e, where G is the normal-state conduc- 
tance of the junction. On the other hand, the minigap 
can be roughly estimated as E^ ^ min(A, -Edwcii), where 
A is the superconducting gap, and E'dwcii = ^idwcU is 
the energy associated with the typical dwell time idwcii 
of an electron in the normal region. Since mesoscopic 
fluctuations for quantum dots and quasi-one-dimensional 
systems are usually small in the parameter Gq/G, one 
can estimate the magnitude of mesoscopic fluctuations 
as Sic ~ eE^/h. For a long diffusive wire, E'dweii ~ Et, 
indicating that results of Ref. 5 are qualitatively correct. 

The scattering matrix approach^^ has proven to be a 
powerful tool for studying coherent electron transport 
in mesoscopic conductors. In the framework of this ap- 
proach, an arbitrary scatterer can be described by a set 
of transparencies in each conduction channel. The main 
transport characteristics of a mesoscopic conductor, such 
as the conductance, shot noise power, conductance of 
the NS junction, can be written as a linear statistics 
on the transmission eigenvalues, J2n fC^n)-^^ Then, the 
corresponding weak localization corrections and meso- 
scopic fluctuations can be determined from the knowl- 
edge of the average density and the correlation function 
of transmission eigenvalues. For a diffusive wire, they 
were obtained in Refs. 16,17. However, the Josephson 
current is generally not a linear statistics on T„'s. It 
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is a linear statistics only for a short (A <^ Et) wire:-*^^ 
/(X) = {eA/2h) ^„r„sinx/[l-r„sin2(x/2)]i/2, where 
X is the superconducting phase difference. In this limit 
one finds for the fluctuations of the critical current at 
zero temperature: ^^'^^ Sic — 0.30(eA/ft). The crossover 
between short and long wires was also considered numer- 
ically within the tight-binding model. 

Recently, the problem of mesoscopic fluctuations of the 
supercurrent in short junctions with weakly transparent 
NS interfaces was considered by Micklitz^" who used the 
machinery of the supersymmetric nonlinear cr-model^^. 
Within that approach, the average Josephson current is 
obtained at the level of the saddle point which corre- 
sponds to the quasiclassical Usadel equations^^, while 
mesoscopic fluctuations can be obtained from fluctua- 
tions around the saddle point. We believe that the func- 
tional (T-modcl approach is a proper tool for studying 
mesoscopic fluctuations of the Josephson current, espe- 
cially in the limit of long junctions where the multiple 
scattering theory fails. 

A couple of experiments with the use of a gated 
semiconductor instead of a normal metallic part have 
been reported. ^'^'^^ Mesoscopic fluctuations were ob- 
served with varying the gate voltage. In Ref. 23, meso- 
scopic fluctuations were shown to follow precisely those 
of G, the latter being estimated from I{V) curves at a 
bias voltage larger than the superconducting gap. On the 
other hand. Sic was systematically found to be smaller 
than the theoretical estimate from Ref. 5. A similar be- 
havior was reported in Ref. 24. 

In this work, we present the derivation of mesoscopic 
fluctuations, SI{x) = [{I'^ix)) - {I{x)fY^^, and weak lo- 
calization correction, Awl^(x), to the Josephson current 
in diffusive SNS junctions within the replicated version 
of the nonlinear cr-model. We ignore interaction effects 
in the normal metal and assume zero magnetic field. By 
expanding the action of the cr-model around its saddle 
point, corresponding to an inhomogeneous solution of the 
Usadel equations, we present mesoscopic fluctuations in 
terms of the soft modes, analogous to Cooperons and dif- 
fusons in the normal state^°'^^. This approach allows us 
to follow the crossover between the regimes of short and 
long wires. 

In general, we verify that for the Josephson current 
through a chaotic dot and quasi-one-dimensional wire, 
SIc/Ic ^ —^whlc/Ic ~ Gq/G, where G is the normal- 
state conductance of the system, and determine the exact 
coefficients in these relations as functions of the junction 
length and temperature. These coefficients are generally 
of the order of 1 , but in some cases (two tunnel barrier 
structures with low barrier transparency) are addition- 
ally suppressed. 

In particular, we find that the approach of Ref. 5 which 
does not take the proximity effect into account systemat- 
ically underestimates rms Ic by the factor of order 3. In 
the limit of long {Et ^ A) quasi-one-dimensional junc- 



tions at zero temperature, we obtain 

We =1.49^, (1) 

which is 2.5 times larger than the prediction of Ref. 5. 
Similarly, for wide {W^.Wy > L) and long {Et < A) 
three-dimensional junctions made of a metallic paral- 
lelepiped of size L X Wy x M^z, we find at T = 0: 

«c = 2.0— (2) 

which is 2.8 times larger than the corresponding result of 
Ref. 5. 

The paper is organized as follows. In Section II, we 
consider the case of a chaotic diffusive dot coupled to 
the superconducting leads through tunnel barriers. This 
model allows us to introduce the method and is simple 
enough to be solved analytically. In Section III, we con- 
sider the case of a diffusive wire with transparent NS 
interfaces. Mesoscopic fluctuations of the critical current 
in two-dimensional (2D) and three-dimensional (3D) ge- 
ometries are calculated in Sec. IV. We discuss the results 
in Conclusion. Technical details are delegated to several 
Appendices. 

II. SUPERCONDUCTOR-CHAOTIC 
DOT-SUPERCONDUCTOR JUNCTIONS 

As a warm-up, in this section we consider a Joseph- 
son junction formed of a "zero-dimensional" chaotic dot 
in contact with two superconducting reservoirs through 
tunnel junctions with conductances Gl and Gr- A pos- 
sible realization of this system would be a short diffusive 
wire, with the Thouless energy Et much exceeding the 
superconducting gap A in the superconductors and with 
the intra-dot conductance Gjv = IttGqEt / 5 (where 5 is 
the mean level spacing in the dot) much exceeding Gl 
and Gfl. As a consequence, the conductance of the struc- 
ture in its normal state is determined solely by the tunnel 
barriers: G = GlGr/{Gl + Gr). 

The average Josephson current in such a system has 
been studied by Aslamazov, Larkin, and Ovchinnikov^^ 
within the tunnel Hamiltonian approach, by Kupriyanov 
and Lukichev^^ with the help of the quasiclassical Us- 
adel equations, and by Brouwer and Beenakker^^ using 
the scattering approach^^ and Random Matrix Theory 
(RMT)^'^. It was found that the amplitude of the super- 
current is controlled by the ratio between the supercon- 
ducting gap A and the "escape" energy Eg ~ G6/Gq. 
The later energy scale is associated to the broadening of 
levels in the dot due to coupling with the leads, play- 
ing here the role of the dwell energy iSdwcii defined in 
Introduction. 

We shall rederive these results in the fermionic replica 
(T-modcl language and then use this formalism to study 
mesoscopic fluctuations of the Josephson current for ar- 
bitrary ratio A/ Eg. A similar approach was very recently 
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followed by Micklitz^" who considered the effect of bar- 
rier transparencies on the average supercurrent and its 
fluctuation in the regime A Eg, within the framework 
of the supersymmetric cr-model. 

In Sec. II A, wc introduce the replica cr-modcl for this 
system. In Sec. II B, we analyze its saddle point solu- 
tion and rederivc the quasiclassical result for the Joseph- 
son current. The fluctuation determinant is calculated 
in Sec. II C. It contains both weak localization correc- 
tions to the supercurrent and its mesoscopic fluctuations, 
which are analyzed in Sees. II D and HE. The results are 
summarized in Sec. II F. 



A. Replica cr-model for a chaotic dot 

The equilibrium supercurrent which flows in a Joseph- 
son junction can be obtained from the free energy = 
—kTlnZ of the system at temperature T: 



2e d 
h dx 



Hx). 



(3) 



where x is the superconducting phase difference between 
the leads. Disorder-averaging is performed in a standard 
way using the replica trick, 



{Hx)) 



-kT lim 

n-+0 r 



1 



(4) 



where {Z'^) can be evaluated as a functional integral 
within the fermionic replica cr-model:^^^"^^ 



VQe 



'S[Q] 



(5) 



The nonlinear cr-model is a field theory formulated in 
terms of the matrix field Q acting in the direct product of 
the replica space of dimension n, infinite Matsubara en- 
ergy space, two-dimensional Gorkov-Nambu space (Pauli 
matrices r^), and two-dimensional spin space (Pauli ma- 
trices (Ti). The Q matrix is subject to the nonlinear con- 
straint = 1 and obeys the charge conjugation symme- 
try: 



Q = Q ^ T1C72Q <J2Tl, 



(6) 



where Q-^ stands for the full matrix transposition. The 
condition (6) is related to simultaneous introduction of 
the Gorkov-Nambu and spin spaces, which renders the 
vectors = {ip^ ,ip*,ipi, —ip*)'^ and be linearly depen- 
dent. The functional integration in Q is performed over 
an appropriate real submanifold of the complex manifold 
defined by the constraints = 1 and Q ~ Q. 

The action of the a-model for a chaotic dot coupled 
to the superconducting terminals via tunnel junctions is 
given by'^^ 



S[Q] = -^t4eTsQ 



E 

i=L,R 



G,6 



Q^Q]■ (7) 



Here, 6 = {i/V)~^ is the mean level spacing in the dot {V 
is the dot's volume and v is the single-particle density of 
states at the Fermi energy per one spin projection), e is 
the fermionic Matsubara energy, and the trace is taken 
over all spaces of the Q matrix. 

In the superconducting reservoirs with the order pa- 
rameters Ae~^^^'^ (left) and Ae'^/^ (right), the matrices 
Qi {i = L, R) are unit matrices in the replica and spin 
spaces, being diagonal in the energy space with the ma- 
trix elements: 

/ X X\ 
Ql,r = [Ti cos - ± T2 sin -j sin + T3 cos 9s, (8) 

where cosOg = cos9s{e) = e/v? + A^. 

Equation (8) is often referred to as "the rigid bound- 
ary condition" . It corresponds to neglecting the inverse 
proximity effect as well as depairing effect in the leads 
due to a finite current density (see, e.g., review 14). 



B. Saddle point: average Josephson current 

We start the analysis of the cr-model (7) with the 
saddle-point approximation, which amounts to neglect- 
ing mesoscopic fluctuations and weak-localization correc- 
tions. The matrix Qo which extrcmizcs the action solves 
the saddle-point equation: 



= 0. 



(9) 



Equation (9) is nothing but the Usadel equation for the 
quasiclassical Green's function for a chaotic dot supplied 
by the Kupriyanov-Lukichev boundary conditions'^ at 
the tunnel interface. This equation can also be obtained 
with the help of Nazarov's "circuit theory" .^^ 

The solution of Eq. (9) proportional to the unit matrix 
in the replica and spin spaces can be easily found: 



where 



]o ~ [ti cos — T2 sin (/)) sin + cos ( 

tan^(.)^ /"^^(^^ , 
tan0 = — — — - tanx, 



(10) 

(11a) 
(lib) 



and 



Esix) 



GI+GI + 2GlGrcosx, (12) 



with Eg = Eg{0). The pole of Qq located at imaginary 
e is related to the minigap i?*(x) in the density of states 
of the normal island.^ In the limiting cases. 



^*^^)-\(A/i^,)i^,(x), 



A<^E, 



9' 



(13) 



9' 



4 



0.8 
0.6 
< 0.4 

* 

0.2 




0.01 



FIG. 1: The minigap vs. A/ Eg in a superconductor- 
quantum dot-superconductor Josephson junction with sym- 
metric tunnel barriers. The two curves correspond to in 
units of A (solid line) and in units of Eg (dashed line). 



while in the intermediate region, A Eg, the depen- 
dence of i?* on X is more complicated. In what fol- 
lows, we will denote i?* = E^{0). Roughly speaking, 
E^: « min(A, Eg), see Fig. 1. 

The action at the saddle point is given by nSo, where 
n is the number of replicas and 



Soix) 



(14) 



with the summation over the fermion Matsubara energies 



Tr{2p+l)kT, and 



Egr 



A^Egixy 



(15) 



In the leading order in G/Gq ^ 1, one can neglect weak 
localization and mesoscopic fluctuations effects (they will 
be studied in the next subsections). In this approxima- 
tion, equivalent to the standard quasiclassical analysis, 
the Gaussian integral near the saddle point yields unity, 
and the average Josephson current can now be obtained 
from Eqs. (3)-(5) as (/(x)>o = i2ekT/h)dS„{x)/dx, 
yielding 



a(x))o 



TkT 



A^ sinx 



(e2 + A2)c.,(x)- 



(16) 



The result (16) is certainly not new. It had been ob- 
tained previously by a number of authors^"'^^"^^. Here 
we simply rederivc this result for completeness. 



Symmetric junction 

The general expression (16) simplifies for symmetric 
barriers: Gl = Gr = 2G. In this case, Eg{x) = 
Eg |cos(x/2)| and Eg = GSIttGq. 

At zero temperature (more precisely, kT E^), the av- 
erage Josephson current is controlled by the ratio A/ Eg. 
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FIG. 2: Quasiclassical results for a superconductor-quantum 
dot-superconductor Josephson junction with symmetric tun- 
nel barriers: (a) the critical phase Xc (in units of vr/2) vs. 
A/ Eg at zero temperature; (b) the critical current Ic (solid 
line: in units of GA/e, dashed line: in units of GEg/e, dotted 
line: in units of GE^/e) vs. A/ Eg at zero temperature. 



When A <^ Eg, the Josephson relation is not sinu- 
soidal: 



(-^(x)>o = — sinxK (sin^ 
e \ 2 



(17) 



where K(x) — F{tt/2,x) is the full elliptic integral of 
the first kind defined as in Ref. 34. The critical current 
Ic — 1.92G'A/e is achieved at a phase difference Xc — 
1.18(7r/2). 

In the opposite limit, A^ Eg, the Josephson relation 
is close to sinusoidal: 



(/(X)>c 



GE„ 



sin X In 



2A 



(18) 



with the critical current Ic ~ (G£^g/e) ln(2A/£'g) at a 
phase difference Xc — 7r/2. 

The crossover for the critical current and critical phase 
at zero temperature and arbitrary relation between A 
and Eg is illustrated in Fig. 2. 

At temperature close to the critical temperature of the 
leads, Tc, the superconducting gap vanishes as A(T) oc 
k[Tc[Tc - T)]i/2. Then for relevant energies e>kT^ 
A{T), Eq. (16) can be simplified, yielding a sinusoidal 
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Josephson relation with the eritieal current 



UnGAyAekT,), A « 

[7C(3)G£;gAV4e^2py2^ A > £;g, 



(19) 



where C(3) w 1.202 is the Rieniann zeta function. 

At intermediate temperatures, Eg <g; kT <^ kTc, we 



find Ic = iGEg/e)ln{2-fA/TTkT), where 7 = e"^ = 
1.781 ... is the Euler constant, in agreement with 
Ref. 25-27. 



C. Gaussian fluctuations near the saddle point 

In this subsection we take Gaussian fluctuations near 
the saddle point into account. Since we will be interested 
in mesoscopic fluctuations of the Josephson current (see 
Sec. HE), we will have to consider the average (Z'^^Z^^) 
of two partition functions calculated at different super- 
conducting phases xi a-nd X2- This average can also be 
expressed in terms of the cr-model: 



{Z 



= / VQe 



-S[Q] 



(20) 



where the only difference with the cr-model described in 
Sec. II A is that now Q becomes an (711 +712) x {ni +^^2) 
matrix in the replica space. Correspondingly, the super- 
conducting Q-matrices in the terminals should be modi- 
fled. Now they are diagonal in the replica space, having 
the superconducting phase difference xi (resp. X2) the 
ni first (resp. 712 last) replicas. With these modification, 
the action of the cr-model has the same form (7). 

At the saddle-point, the matrix Qo extremizing the 
action is diagonal in the energy and replica spaces. Its 
diagonal elements are given by Eqs. (10) and (11), where 
the phase x is set to xi (resp. X2) if 1 ^ a ^ '^1 (resp. 
ni < a < ni + , where a is a replica index. 

In order to study fluctuations near this saddle point 
we write matrices close to Qq as 



= t/U(l + W + +...)U. 



(21) 



Here, the matrices A and U should be chosen in such 
a way that in the absence of fluctuations, at Vt^ = 0, 
Eq. (21) reduces to WAU = Qq. As usually, one has to 
require 

{A,W} = 0, (22a) 

and impose the constraint following from Eq. (6): 

W = -~W, (22b) 

and the requirement of convergency of the cr-model on 
the perturbative level: 



W 



-W. 



(22c) 



The form of the parametrization (21) is standard, while 
the choice of the matrix A is a matter of convenience. 



A possible choice could be the metallic saddle point, ■^^ 
Am = T3 sign(e) [the hmit of Eq. (8) at e A]. Here we 
adopt an alternative choice proposed by Ostrovsky and 
Feigel'man'^^: 



A: 



n, 



(23) 



corresponding to the superconducting saddle point (8) at 
zero energy and % = 0. With this choice of A, the unitary 
matrix U in Eq. (21) is given by 



U 



g-ir2(7r/4-e/2)g-iT3(</./2)^ 



(24) 



The choice of the parametrization with A = ri has two 
technical advantages: (i) the solution of the constraint 
(22a) is independent on energy, and (ii) the constraint 
(22b) can be easily resolved. 

A general parametrization of the matrix W satisfying 
the constraint (22a) is given by: 



W„ 



T2C„ 



(25) 



where m = (e, a) and n — (e', b) encode both the energy 
and replica indices, and dmn and Cmn are 2x2 matrices 
in the spin space which can be expanded in the Pauli 
matrices as 



d ~ do + dcr, 



Co 



(26) 



The variables do and cq (d and c) will be referred to 
as singlet (triplet) modes. They play the same role as 
diffusons and Cooperons in a normal metal, describing 
soft diffusive excitations on top of an inhomogeneous 
proximity-induced state. Note that contrary to diffusive 
modes in a normal metal, these d- and c-modes are gen- 
erally coupled in the presence of a supercurrent in the 
normal region, cf. Sec. III. 

Equations (22b) and (22c) yield 



do = 
Co 



T t 

c = c = c' . 



(27a) 
(27b) 



Here, transposition acts in the replica and energy spaces. 
In terms of the matrix elements, Eqs. (27a), e.g., read: 

(c^o):,'^, = (do)'J'.^-e = (28a) 
[d^rA, = -{d,t^,,,_, = Hd:)%- (28b) 

Independent integration variables for the singlet rfo mode 
can be chosen, e.g., as: 



(rfo) 



ah 



if a > 6, e > 0, 



(do):,":, eC, ife>|e'|>0, 
(do)r-. € C, if e > 0, 
{do)Z G if e > 0, 

and analogously for the triplet Ci=i.2,3. Independent in- 
tegration variables for the triplet di=i,2.3 modes can be 
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chosen, e.g., as: 



aa 



if a > 6, e > 0, 
if e > |e'| > 0, 
if e > 0, 



and analogously for the singlet cg. 

Expanding the action in the Gaussian approximation 
over fluctuations near the saddle point, one finds in gen- 
eral: 



mn i—0 



(29) 



where Amn is a symmetric [it can be symmetrized using 
relations (27)] matrix in the (d, c)-space with the simple 
block structure in the replica space: 



At 



AX1X2 
j\X2Xl 

j^X2X2 



if a, 6 < TT-i, 
if a < TT-i < 6, 
if 6 < Til < a, 
if ni < a, b. 



(30) 



The matrix A does not depend on the spin index i since 
the spin in conserved. Such a dependence will arise if one 
takes magnetic impurities or spin-orbit interaction in the 
normal region into account. 

Due to the absence of the {di)'^''_^ and (co)""!^ modes, 
the matrix A'^°_^ should be diagonal in the (d, c)-space: 



XX 











(31) 



We will sec below that this is indeed the case [cf. Eqs. (35) 
and (62)]. 

Integration over independent variables of the d and c 
modes gives the fluctuation determinant: 



{Z"^Z"'^) — e^"^^oixi)-n2Soix2) ^ 
where Sq{x) contains the WL correction: 



dd 



(32) 



(33) 



[we write this expression in the most general way assum- 
ing that (A^^^)'^'^ and {A^^^y^ might be operators, as in 
Sec. Ill; for a chaotic dot considered in this section, the 
trace in Eq. (33) can be omitted], whereas the prefactor 
Ai accounts for mesoscopic fluctuations: 



1 



(det^^,'-,'^^)"^ (detA^^^,'^^)'"^"^ (detA^,^,'^^)"= 

(34) 

In Eq. (34), the product over e and e' should be taken over 
all Matsubara energies and we have omitted the factors 



which are equal to 1 in the replica limit ni.2 — > and do 
not depend on Xi,2- 

For a superconductor-quantum dot-superconductor 
junction considered in this Section, the matrix A can be 
easily found by expanding the action (7) in W with the 
help of Eqs. (21)-(25): 



A3 



(35) 



where coeix) is defined in Eq. (15), and Eq is the identity 
matrix in the (d, c)-spacc. 

The weak localization correction to the Josephson cur- 
rent and its mesoscopic fluctuations are discussed in the 
next Subsections. 



D. Weak localization correction 

The weak localization correction to the Josephson cur- 
rent, AwL-^(x) = (^(x)) — (^(x))o: can be found with the 
help of Eqs. (3), (4), (32) and (33): 



Awl/(x) 



ekT a^^^^^(Af_J 



dd 



h dx 



{A 



XX \ 
e. — e) 



(36) 



Since for a superconductor-quantum dot- 
superconductor junction, the matrix A [see Eq. (35)] 
acts as a unit matrix in the (d, c)-space, there is no weak 
localization correction to the Josephson current. 

Note that there is no weak localization correction to 

the conductance of a normal metal-quantum dot-normal 

metal junction in the lowest order in the tunnel barriers 
eitheri5'36^ 



E. Mesoscopic fluctuations 

With the help of the replica trick and relation (3), the 
current-current correlation function can be expressed as 



{iixi)nx2)) 



2ekTV (Zl^lZl^-) 
lim 



dxidx2 ni,n2^0 nin2 



(37) 

Making use of Eqs. (32) and (34), we get the gen- 
eral expression for the cumulant ((/(xi)/(x2))) ^ 
(/(Xi)/(X2))-(/(xi))(/(X2)): 

((/(x.)/(X2)))^-8(f:)^^|:trlnA-. 

(38) 

Since A^^,^^ given by Eq. (35) appears to be diagonal 
in the (d, c)-space, evaluation of Eq. (38) for the quantum 
dot geometry is trivial and we get: 



(a(xi)/(X2))) 



AekT 



E 



1 [^£(Xl) + t^e'(X2)]^ 9X1 l9X2 



(39) 
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FIG. 3: (a) Mesoscopic fluctuations of the critical current, 
Sic = {vaxlc)^^^ (solid line: in units of eA/h, dashed line: in 
units of eEg/h, dotted line: in units of eEt^/K) vs. A/Eg at 
zero temperature, (b) The ratio GSIc/GqIc vs. A/Eg at zero 
temperature. 



In the limits T = and A Eg, Eq. (39) had been 
recently derived by Micklitz^" using the supersymmetric 
cr-model approach. 

As shown in Ref. 37, mesoscopic fluctuations of the 
critical current, 5Ic, can be obtained as Sic ~ 61{xc)- 



Symmetric junction 

We consider now symmetric junctions. 
At zero temperature, wc find: 

^ f0.396(eA/ft), if A « 
' [eEg/nh, if Eg <^ A, 

while the result for an arbitrary ratio A/ Eg is plotted in 
Fig. 3. 

A peculiar feature of the T = mesoscopic fluctuations 
given by Eq. (39) is a finite limit at x ~* ""^ 

nfi A + Eg h 

while the Josephson current must vanish exactly at x 
.^ 20,37 'pjjjg ig related to the breakdown of the Gaus- 




FIG. 4: Mesoscopic fluctuations of the Josephson current, 
SI{x) vs. phase x at zero temperature (solid line) and at tem- 
perature kT — E*/10 (dashed line). The curves are plotted 
for different ratios A/Eg-. (top), 0.2, 0.5, 1, 2, 5 (bottom). 



sian treatment of fluctuations at % tt. A more care- 
ful analysis of the general nonlinear action (7) should 
restore the exact relation /(tt) = 0, resulting in van- 
ishing fluctuations for x ^ At a small tempera- 
ture, 6I{'!t) = 0, and the Josephson current decreases 
from its typical value (41) to in the small phase range 
|x-7r| < kT/E^, see Fig. 4. 

Close to the critical temperature, Tc, when A(r) — > 0, 
we find 



SI 



2ekT 



E 



A'iT)El 



e,e'>0 



e^e'^ie + e' + 2Egr{e + Eg)ie' + Eg) 



(42) 



Thus, 



(eA2(r)/8ftfcT,), 



if kTc^Eg, 



' \O.OlOeA'^{T)E^g/hk^Tj^, ii Eg <^ kT^. 

At intermediate temperatures such that Eg <^ kT <JC 
kTr, we find 



SI 



2ekT , 



1 



^ eeUe + e') 

e.c'>0 ^ ' 



l\2 ■ 



(44) 



Thus, Sic ^ 0.l2{eE^/hkT). This result was also ob- 
tained by means of the fourth order perturbation the- 
ory in the tunnel Hamiltonian connecting the wire to 
the leads, plus random matrix theory on the statistics of 
eigenstates in the lead.'^^ 



F. Summary of results 

The results of our study are summarized in Table I. 
We note that, quite generally, the critical current at low 
temperatures is related to the minigap in the dot at x = 



8 



TABLE I: Summary of results for a quantum dot contacted symmetrically to superconducting leads {Gl = Gr — 2G). 



elc/G 



Jix) 



hSIc/e 



gSIc/Ic 



A^Eg,T = 
£9 < A, T = 



1.18 
1 



arbitrary Eg/ A, T = Fig. 2(a) 
kT < kT, < Eg 1 
Eg<.kT< kT, 1 
E„«.kT^ kTc 1 



^ g2^[26,271 

Eg\n{2A/E, 



\[271 



Fig. 2(b) 
7rA2(r)/4A:rJ2'^l 
Q.21?,EgA^{T)/{kTcf f^*^' 
£;<,In(27A/7rfcr)P5-2^1 



0.396A'^''1 

Eg/^ 

Fig. 3(a) 
A2(r)/8A;Tc 

o.oioA2(r)£;^/(fcT,)3 

Q.12El/kT 



0.648 
l/ln(2A/iJ,) 
Fig. 3(b) 
1/2 

Q.lbEg/kT^ 

0.38^9 / [kT ln(27A /nkT)] 



Here, g = G/Gg, A(T) = [(87r2/7C(3))fc2Tc(Tc - T)]i/2 and 7 = 1.781 . . 



0, E^, ^ min(£'g,A), through an Ambcgaokar-BaratofF 
like formula: Ic ~ GE^je^ where G is the normal-state 
conductance of the system (up to a logarithmic factor for 
E, « A). 

In a normal metallic double-barrier structure, conduc- 
tance fluctuations are "universal": &G = Gq/2^ while 
weak localization corrections vanish at infinitely small 
mean level spacing in the dot^^'^^. 

In an S-QD-S junction, we find that the WL correc- 
tion to the Josephson current also vanishes, while the 
amplitude of mesoscopic fluctuations can be estimated 
as 5Ic/Ic ~ 5G/G, provided kTc Eg. In the opposite 
limit, for poorly conducting barriers, mesoscopic fluctu- 
ations are additionally suppressed. 

Note that in the limit kT < kTc *C Eg, we have 
SIcJIc = SG/G. This should be attributed to the fact 
that in this limit the Josephson relation becomes sinu- 
soidal and proportional to the conductance of the system. 
Indeed, using Eq. (15) from Ref. 18 with A(T) < kT, we 
find /(x) = ^cSinx, where the critical current is propor- 
tional to the exact conductance before disorder averag- 
ing: /e = (eA2(r)/4;ifcr)E„T„ = nA^iT)G/4ekT. 



III. SUPERCONDUCTOR-NORMAL 
METAL-SUPERCONDUCTOR JUNCTIONS 

We turn now to the case of a diffusive metallic wire 
connected to superconducting leads by transparent in- 
terfaces (the interface resistances are much lower than 
the resistance of the wire in the normal state). 

Proximity effect in such a geometry can be described 
by the diffusive replica cr-model with the action^^'^^ 

S[Q] = ^ y drtr [D{VQf - AersQ] , (45) 

where D is the diffusion coefficient. At the boundaries 
with superconductors, we require that Q matches with 
the Qi matrices in the leads given by Eq. (8). 

In this Section we concentrate on a quasi-one- 
dimensional geometry, when the length of the wire, L, 
much exceeds its transverse dimensions. Then, the spa- 
tial dependence of Q is reduced to the dependence on the 
coordinate x along the wire, which will be measured in 



units of L. The action (45) can be written as 

^[Q] = f'^ [(^Q)' - 4eT3Q] , (46) 

16Gq 7_i/2 

where Gm = '2'kGqEt / 5 is the normal-state conductance 
of the wire, and e = c/Et stands for Matsubara energies 
measured in units of the Thouless energy, Et = HD/L^. 

Following the same line as in Sec. II, we derive the 
quasiclassical Josephson relation as the saddle point of 
the action (46) in Sec. Ill A. We find that the amplitude 
of the critical current is controlled by the ratio between 
A and Et- Then, we express mesoscopic fluctuations 
and WL correction to the Josephson current in terms of 
Gaussian fluctuations in the vicinity of the non-uniform 
saddle point in Sec. IIIB. Results are summarized in 
Sec. inc. 



A. Josephson relation 

At the saddle point, the matrix Qa{x) which extremizes 
the action (46) solves the Usadel equation: 

-V(QoVQo)+e[r3,Qo] =0. (47) 

With the parametrization (10), this equation can be re- 
duced to two differential equations: 

V(sin2 6lV(/)) = 0, (48a) 
V^e-2e sin 6 - {Vcj)f sin 6I cos 6I = 0, (48b) 

with the boundary conditions d{±l/2) — 9s and 
(/)(±l/2) = ±x/2. The Usadel angles at energies ±e are 
related by the symmetries: 

= 7r-0e, 0-. -0.. (49) 

According to Eq. (48a), J — siv? 9\/(j) (usually referred 
to as the spectral current) is constant along the wire. 
Integrating then Eq. (48b), we obtain: 

(yef - 4e (cos 0(0) - cose) + - -4-. (50) 

sur t'(O) sm 

This last equation can be solved in quadratures, the im- 
plicit solution being given in terms of elliptic integrals. 
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Its form is rather cumbersome and can be found in Ap- 
pendix A. In a general case, J and ^(0) at given boundary 
conditions can be found only numerically. Note however 
that their determination requires solution of only an al- 
gebraic rather than differential equation. 

Equation (50) can be used to determine the minigap 
induced in the normal region due to the proximity effect. 
To this end, one has to perform analytic continuation 
e — > ~iE, and search for the energy when cos 9 first 
acquires a finite real part, indicating a nonzero density of 
states above . Determination of the minigap simplifies 
for X = 0. To solve the Usadel equation we introduce 
then 9 = 7t/2 + and get the relation between E and 
the value ^^{O) at the center of the wire: 



^(0) 



dip 



ah{E/A) ^sinh'i/;(0) — sinh?/' 



(51) 



Equation (51) establishes a relation between E and "0(0). 
The value of E^, should be determined from the condition 
that Eq. (51) ceases to have real solutions for ip{0). For 
short wires (A <^ Et), iJ* « A, whereas for long wires 
(A » Et), E^ w 3.122£'t,^^'''° in accordance with the 
general relation ^ min(£'T,A). The dependence of 
i?* on A/Et is shown in Fig. 5. 

The action at the saddle point Qq{x) is given by nSo, 
where 



So 



G 



N 



4Gc 



^ j dx [{V9f + [Vpf sin^ 6* - 4£ 



Taking the derivative with respect to Xi integrating by 
parts, and using the boundary conditions at x — ±1/2 we 
get dy^So = {Gn/'2.Gq)Y^^J, yielding the quasiclassical 
expression for the supercurrent: 



(/(x))o=7rfcr^^J. 



(53) 



Below, we evaluate the current (53) in different cases. 



The critical current for a short diffusive wire {Et ^ A) 
has been known for a long time.^^ In this case, the term 
proportional to e in Eq. (48b) can be neglected, and the 
Usadel equation can be solved exactly: 



9{x) 



arccos 



cos 0(0) cos — 



Jx 



sm9{Q) 



where 



sin 61(0) = 



sm f , cos ' 



: sin^ f 



J = 2 sin 6'(0) arcsin I sin 9s sin 



(54) 



(55a) 



(55b) 



Calculating the supercurrent with the help of Eq. (53) at 
T = we get" 

/(x) '^GnA ^ arctanh ^sin . (56) 

The Josephson relation is nonsinusoidal, see Fig. 6. The 
critical current Ic = 2.082GArA/e is achieved at the crit- 
ical phase Xc — 1.255(7r/2). 



2,. Arbitrary wire at temperature close to Tc 

At temperatures close to the critical temperature in 
the leads, Tc, the superconducting gap A in the leads 
is smalh A(r) = [{8tt^ /7C{3))k^Tc{Tc - T)]^^^ < kT. 
Thus, at relevant energies e ^ kT 3> A, the anomalous 
(Gor'kov) component of the Green function is small, and 
the Usadel equation (48b) can be linearized with respect 
to sm9. Its solution has the form 



sin6'(a;)e 



A sinhK(i — x) 

g ' 2 i 

lei sinhK 



A ,x sinhK(i + 

H — — 

|e| sinhft 



x) 



(57) 
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FIG. 8: Josephson relation I{x) (in units of GnEt/s) for a 
long wire at zero temperature. 

considered in Sec. IIIAl. For arbitrary A/Et, the crit- 
ical current and critical phase obtained numerically are 
shown in Fig. 7.^'^ 

Specifically, for long junctions {Et ^ A) the Joseph- 
son relation is still highly nonsinusoidal (see Fig. 8), with 
the critical current 1^ = 10.83GnEt /e achieved at the 
critical phase Xc = 1.271(7r/2).i3 



FIG. 7: (a) The critical phase Xc (in units of tt/2) vs. A/Et 
at zero temperature, (b) The critical current 7c (solid line: in 
units of GivA/e, dashed line: in units of tiGnEt /e, dotted 
line: in units of Gjvi5*/e) vs. A/Et at zero temperature. 



where k = y^jej = \j2\e\/ET- Calculating the spectral 
current with the help of Eqs. (Dl) and (D2), we get 



J 



sinh / 



■ smx- 



(58) 



Close to Tc, junctions are classified as short or long 
depending on the ratio between kTc and Et- For short 
junctions [A(r) <C kT < kTc ^ Et], the supercurrent is 
given by: 



i{x) 



4e kf 



smx- 



(59) 



For long junctions [A{T),Et kT < kT^], the super- 
current is exponentially suppressed: 



Hx) 



SG 



N 



l2iTkT 



: exp 



sinx. (60) 



e ^2TikTET \ V Et 
In both cases the Josephson relation is sinusoidal. 

3. Arbitrary wire at zero temperature 

At T = the Josephson relation /(x) depends only 
on the ratio between A and Et- Short junctions were 



B. Mesoscopic fluctuations and weak localization 
correction 

Fluctuations near the saddle point Qo{x) can be 
parametrized using Eqs. (21), (24) and (25). Substitut- 
ing these expressions into the action (46) and expanding 
to the second order in the modes d{x) and c{x), we get 
similar to Eq. (29): 



mn 2—0 

"■(61) 

The form of the operator A{x) is quite cumbersome, see 
Eq. (B2) in Appendix B. Apart from the second deriva- 
tive with respect to x it contains also the first derivative. 
The latter can be eliminated by a proper unitary trans- 
formation [Eq. (B5)] mixing the d and c components of 
fiuctuations. After such a rotation the matrix A acquires 
the form"*^ 

Aran{x) = -V^ + Uran + Pnin COs(ry,„ + r/„)I]3 

4- pran sin(7?™ + r/„)I]i, (62) 

where 



ocmn = Em COS 6*™ - -[(V6'„)2 -|- (sin 6',„ V0,„)^] 

+ e„cos0„ - i[(V0„)2 ^ (sin0„V0„)2], (63) 



JV(V&^)' + (sm0™V0„)2 



X V(V^n)' + (sin0„V<^„)2, (64) 
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and the odd function rjm {x) can be obtained by integrat- 
ing the relation 



(V0„O2 + (sin0™V0™)2 



(65) 



Once the operator A is known, one can use the gen- 
eral Eqs. (36) and (38) to calculate the weak localization 
correction and mesoscopic fluctuations of the Joscphson 
current. The only difference compared to the S-QD-S 
case is that now A{x) is an operator in the real space, 
and the determinant should be calculated with respect 
to spacial coordinates as well. In a general situation, 
for arbitrary and temperatures, this can be done 

only numerically, following the procedure outlined in Ap- 
pendix A. 

Several cases where a simple form of the operator 
A{x) allows for analytic solution are discussed below. In 
Sec. IIIB 3, wc present the results of numeric solution for 
an arbitrary A/Sy at zero temperature. 



1. Short wire 

We start with the simplest case of a short wire, A <C 
Et- This limit was analyzed previously in Refs. 16 and 
17 using the scattering-matrix approach to the Josephson 
current. 

For short wires, the term proportional to e in the 
Usadel equation can be neglected. Then, according to 
Eq. (50), {VOra? + (sin0^V</.™)2 = 4C^, where = 
Jm/2sin0m(O) is constant in the wire. With the same 
accuracy, Eq. (65) guaranties that rjmix) = 0. As a re- 
sult, the operator A in Eq. (62) becomes diagonal in the 
(d, c)-space, with its diagonal elements being diffusion 
operators: 



Ar, 



-V2 - (C,n - Cnf 

- (C,„ + C„)2 



Since fluctuations must vanish in the reservoirs, the 
eigenvalues of — are Ti^p^ (p = 1, 2, . . . ), and the de- 



terminant of A^^/^^ involved in Eq. (38) can be readily 
obtained. As a result, the current-current correlation 
function acquires the form 



{{Hx)iix'))) = 



dxdx 



sin(C - C") sin(C + C") 



C-C 



c + c 



(66) 



where C 



3in[(A/VA2 



)sin(x/2)] and C" is 



given by the same expression with e ^ e' and x ~' x' ■ 

At zero temperature, the sums over Matsubara ener- 
gies reduce to integrals, and our expression for SI^{x) = 
((/^(x))) = var/(x) becomes equivalent to the result of 
Refs. 16 and 17. The equivalence is explicitly demon- 
strated in Appendix C. For small Xj we have an expan- 
sion in powers of sin^(x/2): 



siHx) 



eA 
IT 



^ sin^ X 
30 



62 . 2X 
— sm — 

63 2 



(67) 



The whole curve SI{x) at T = is plotted in Fig. 9. 
Mesoscopic fluctuations of the critical current. Sic = 
61{xc), are characterized by Sic = 0.30eA/?i.^^'^^ 

At zero temperature, the magnitude of mesoscopic 
fluctuations remains finite at x = tt, SI{tt) = eA/nh, 
contradicting the general symmetry requirement of van- 
ishing /(x).^°'^^ This is related to the breakdown of 
the Gaussian treatment of fluctuations at x ^ A 
more careful analysis of the general nonlinear action 
(46) should restore the exact relation /(tt) = 0, result- 
ing in vanishing fluctuations for x ^ "'•^'^ The region 
where the exact treatment should modify the result ob- 
tained in the Gaussian approximation can be estimated 
as Ix — ■"'I ^ Gq/Gn <C 1. Finite temperatures ren- 
der SI{tt) ~ 0, and shrink the region of strong non- 
Gaussian fluctuations near x = ""j which disappears at 
A:T» (GQ/G^)|x-_7r|A. 

The weak localization correction can be evaluated with 
the help of Eq. (36). Here, only c-modes contribute to 
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the result, and 



, efcT 9 sin2C 
Awl/(x) = — — Z^ln- 



h dx 

At zero temperature, we obtain: 



2C 



AwL/(x) = -^sinx(l + ^sin^| 



(68) 



(69) 



in agreement with Eq. (3.33) of Ref. 16. The whole 
dependence Awl-^(x) at T = is shown in Fig. 10. 
The WL correction to the critical current is Awl^c = 
-0.266eA/n. 

At zero temperature, the WL correction (68) is discon- 
tinuous at X = ^' AwL-f(7r T 0) = ^eA/Ah, which is an 
artefact of the Gaussian approximation employed. The 
situation here is completely analogous to the situation 
with mesoscopic fluctuations discussed above. 



2. Arbitrary wire at temperatures close to Tc 

In Sec. Ill A 2 we have seen that calculation of the qua- 
siclassical Josephson current simplifies at temperatures 
close to Tc- The same occurs for mesoscopic fluctua- 
tions and WL correction. In order to calculate them, 
we decompose the operator A = Aq + V in Eq. (62) 
into a sum of the diffusion operator in the normal state, 
Aq = -V2 + |ei| + |£2|, and the perturbation F = ©(A^). 
Expanding tr In A in powers of V we rewrite Eqs. (36) and 
(38) as 

Awl/(x) = -^^Y.^tG±sK%, (70) 



h dx 



((/(Xi)/(X2)))=4 



dxidx2 



^Tr(Gf)^ 



(71) 



where G ^ Aq is the Green function of the diffusion 
operator: [-'V'^ + \£i\ + \e2\]Geie2{x,y) = S{x-y), and 
Tr implies tracing over coordinates as well. Explicitly, 



^6162 i^i y) 



sinhA(i +m) sinhA(i - M) 



(72) 



A sinh A 

where m = min(a;,?/), AI = max(a;,?/), and A = 

Vkil + k2|. 

We start with the weak localization correction. Ex- 
panding 



V{x) = F°(x) -t- V^{x)ti + F^(x)E3, 



(73) 



with V'^{x) given by Eqs. (D4)-(D6), we get for the ex- 
pression entering Eq. (70): 



TrGSgV; 



XX 



1/2 

2 / dxGe.-e{x,x){V^{x)) 

'-1/2 



XX 



(74) 
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FIG. 11: (a) Weak localization correction to the critical cur- 
rent — Awl/c (solid line: in units of eA//i, dashed line: in 
units of neEg/h, dotted line: in units of eE^/h) vs. A/Et 
at zero temperature, (b) The ratio —AwlIc/Ic (in units of 
Gq/Gn) vs. A/Et at zero temperature. 



Evaluating the integral using Eqs. (72) and (D6) we ob- 
tain the weak localization correction: 



Awl/(x) 



ekT . ^ A^ cosh2K , , „ , 

■ smx 7, — 3 — ~ tanh2K), 



8h 



sinh'' 



K 



where k = ^/2\e\ = y2|e|7&r. 
For short junctions (fcTc ^ Et), 

„ A2 

Awl/(x) = - 



i2;ifcr 



smx, 



(75) 



(76) 



and Awl/(x)//(x) = ~Gq/3Gn. 
For long junctions {Et ^ kTc), 



Awl/(x) = - 



eA2 



TT h^/2TTkTET 



-^2TrkT/ET 



sinx, (77) 



and Awl/(x)//(x) - -Gq/2Gn- 

Mesoscopic fluctuations are calculated in Appendix D 
with the help of Eq. (71). The result has the form 



((^(Xi)^(X2))) = 4 ( ^ ) sinxi sinx2 



E 



A* / Kl 



e\e\ \sinhKiy \sinhK2 



K2 



(78) 
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FIG. 12: (a) The amplitude of mesoscopic fluctuations of 
the critical current 5Ic (solid line; in units of eA/h, dashed 
line: in units of neET/h, dotted line: in units of eE^/h) vs. 
A/Et at zero temperature, (b) The ratio 5Ic/ Ic (in units of 
Gq/Gn) vs. A/Et at zero temperature. 



where the function Te^e^ defined in Eq. (DIG). 

For short junctions (fcTc ^ -Et), one can take the hniit 
Ki,K2, A — > and get ^eie2 — 1/30, leading to 



{{i{xi)i(a))) 



eA2 



120 \hkT 



(79) 



This result could have been deduced already from the 
exact result (66) for short junctions. The relative fluc- 
tuations of the Josephson current are SI{x)/I{x) = 

y^2/T5 Gq/Gn ■ 

For long junctions {Et ^ kTc), the sums in Eq. (78) 
are dominated by the lowest Matsubara frequencies 
£1,62 = ±7rT, and Ttjt,ttT = e^'*/128K^ with k = ki = 
K2- Hence 



{{IiXi)IiX2))) 



4e2A4 



TT^h?kTET 



-^f,TTkT/Ej 



smxism X2, 



(80) 

and the relative mesoscopic fluctuations are 
«(X)//(X) = Gq/2V2Gw. 



3. Arbitrary wire at zero temperatures 

Finally, we discuss the case of an arbitrary wire at 
zero temperature, when the problem can be solved only 
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FIG. 13: Weak localization correction —Awhlix) (in units of 
eEr/h) vs. x for long wires at zero temperature. 




FIG. 14: Mesoscopic fluctuations SI{x) (in units of eEr/h) 
vs. X for long wires at zero temperature. 



numerically. Here, the amplitude of mesoscopic fluctua- 
tions, 51 (x), and WL correction, AwL-^(x), are functions 
of two parameters: the ratio A/Et and the phase differ- 
ence X- 

The dependencies of the WL correction to the crit- 
ical current and its mesoscopic fluctuations on A/Et 
are presented in Figs. 11 and 12, respectively. One can 
clearly see the crossover from the short to long limits at 
A ^ 5Et- Note that the relative magnitude of meso- 
scopic fluctuations shown in Fig. 12(b) is nearly insen- 
sitive to the wire's length and can be approximated by 
SIc/Ic - O.UGq/Gn. 

For long junctions (Et ^ A), the energy scale for 
mesoscopic fluctuations is set by the Thouless energy: 
6Ic = lAmeEr/h and Awl^c = -1.754e£:T/^- The x- 
dependences of the WL correction and mesoscopic fluc- 
tuations in this regime are shown in Figs. 13 and 14, 
respectively. These plots are inaccurate in the vicinity of 
X = TT, cf. discussion in Sec. IIIB 1. 

Note that our result for the magnitude of mesoscopic 
fluctuations of the critical current for long wires is 2.5 
times larger than the prediction of Ref. 5. 
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TABLE II: Summary of results for a wire with conductance Gjv contacted to superconducting leads. 



Xc/W2) 



hSI^/e 



gNSIc/Ic 



QNAwhlc/ Ic 



A<St, T = 1.255 

£t < A, T = 1.271 

arb. Et/A, T = Fig.7(a) 

kT < kTc < Et 1 
Et <.kT < kTc 1 



2.082A1"' 
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1.49QEt 
Fig. 12(a) 
A^(^)/^/l20fcT 
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0.45 
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Fig. 12(b) 
\/27T5 
1/(2^2) 



-0.266A'^^1 
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Fig. 11(b) 
-A'^{T)/l2kT 
~4A'^{T)/tv£ 



-0.401 
-0.509 
Fig. 11(b) 
-1/3 
-1/2 



Here, gjv = Gn/Gq, A(T) = [{Sw^ /7(;{3))k'^Tc{Tc - T)]i/2^ and £ = ^2-KkTET exp( v'Stt/cT/Et). 



C. Summary of results 

The results of this Section arc summarized in Ta- 
ble II. The critical current at low temperatures is re- 
lated to the minigap in the normal region at x = 0, 
^ min(i?T, A), through an Ambegaokar-BaratofF like 
formula: Ic ~ GnE^/c. We see that, quite gener- 
ally, mesoscopic fluctuations and WL correction to the 
Josephson current arc suppressed by the factor Gq/Gn- 
Sic — AwL^c ~ eE^:/h, where the precise coefficients 
in this expression depend on the wire's length and tem- 
perature. Note that contrary to the S-QD-S cases con- 
sidered in Sec. II, these coefficient arc always of the or- 
der of 1. This should be compared to the behavior of 
the conductance. In a normal metallic wire, conduc- 
tance fluctuations and its weak localization correction 
are also suppressed by the same factor Gq/Gn, but 
the coefficients are "universal": 6Gn = \/2/15 Gq and 
AwlGat = —Gq/3. 

Note that for short wires at T close to T^, the relative 
WL correction and mesoscopic fluctuations of Ic coincide 
with those of the conductance. This is due to the relation 
I{x) oc Gsinx, see discussion in the end of Sec. II F. 



IV. MESOSCOPIC FLUCTUATIONS IN 2D 
AND 3D GEOMETRIES 

A. Mesoscopic fluctuations of the Josephson 
current through a 2D electron gas 

In Ref. 23, mesoscopic fluctuations of the critical cur- 
rent had been measured in junctions formed by a two- 
dimensional electron gas whose width W was much longer 
than the distance L between the superconducting elec- 
trodes. In this case, the results of Sec. Ill cannot be ap- 
plied since transverse diffusive modes become relevant. 
In this Section we take them into account and calculate 
rms 5Ic for the two-dimensional geometry. 

Mesoscopic fluctuations can be obtained with the help 
of the general formula (38), where now the operator A 
reads 

i(2°)(x,2/) = -|^+i(x). (81) 



Here, A{x) stands for the operator (62) for the quasi-lD 
wire, and -l/2<a;<l/2(0<y< W/L) is the dimen- 
sionless coordinate along (perpendicular to) the junction. 
The eigenvalues of the operator (81) are given by 

^mA = — 5 h Ai, (82) 

where 77 = W/L, and Ai are the eigenvalues of the oper- 
ator A{x) (we suppressed the indices X1jX2 and e,s' for 
brevity). For wires, only the zeroth transverse diffusive 
mode with m = was relevant, whereas for the film ge- 
ometry one has to sum over all m's. As a result, wc can 
express the current cumulant in terms of the spectrum of 
the operator A{x): 



mxi)iix2))) 

where 




is a generalization of the function In A relevant for wires 
to arbitrary ratios 77 = W/L. 

We analyze Eqs. (83) and (84) in the experimentally 
relevant limit^'^ of wide {W S> L) and long {Et = 
HD/L'^ ^ A) junctions, at small temperatures {kT <C 
Et)- Then F(A) w r]VX, and the spectral sum should 
be calculated numerically, analogously to the quasi-one- 
dimensional situation, see Sec. Ill B 3. We obtain for the 
rms of the critical current fluctuations (at Xc = 1.277r/2): 

cEt fW 
5Ic^l.5^J^. (85) 
n \ L 

A similar equation with the prefactor 2.2 was used in 
Rcf. 23 [see Eq. (2) there] as the theoretical estimate 
for the critical current fluctuations. Though the two- 
dimensional case was not considered by Altshuler and 
Spivak^, it was claimed that Eq. (2) of Ref. 23 can be 
obtained from the three-dimensional result of Ref. 5. We 
could not follow the derivation of Eq. (2) in Ref. 23 but 
would like to emphasize that even the three-dimensional 
result for 5Ic from Rcf. 5 is overestimated in Eq. (1) of 
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Ref. 23 by the factor of tt^ due a different definition of 
the Thouless energy. So we expect that the 2D result for 
5Ic obtained within the Altshuler-Spivak approach would 
have the form (85) with the prefactor smaller than 1. 

Thus, even though the approach of Ref. 5 generally 
underestimates the magnitude of critical current fluctua- 
tions, the theoretical prediction of Ref. 23 overestimated 
them by the factor of 1.5. 

With the help of Eq. (85) we have a better explanation 
of experimental results (in the weakly localized regime I) 
from Ref. 23: 51^^''°' = 45 nA, 51^"^"' = 15 nA. StiU, 
discrepancy by the factor of 3 remains. 

For completeness, we present also the result for wide 
{W > L) and short (A < Et) junctions at T = 0. It 
can be obtained using the spectrum found in Sec. IIIB 1. 
We get for the rms of the critical current fluctuations: 



eA W 
= 0.26^J-^. 
h \ L 



3D case 



(86) 



Finally, we mention that the above results for the 2D 
case can be easily generalized to the 3D case. To be spe- 
cific, we consider the limit of wide junctions, when both 
transverse dimensions are large: Wy,Wz ^ L. Then, 
mesoscopic fluctuations can be found with the help of 
Eq. (83) with F(A) = {WyWz/4TTL'^)X\nM/X, where 
M ~ L/Z is the high- momentum cutoff [it drops from 
the answer since tr A/dxdx' = 0, see Eq. (62)]. 

For long junctions {Et <C A) at T = 0, numerical 
integration leads to Eq. (2). Note that the result of Ref. 5 
obtained in the same limit contains the numerical factor 
■\/l5C(5)/7r3 = 0.71 instead of 2.0. Again, we see that 
the approach of Ref. 5 underestimates the magnitude of 
mesoscopic fluctuations of Ic by a factor of 2.8. 

For short junctions (A <C Et) at T = 0, we obtain 



a/e = 0.28-^^ 



CONCLUSION 



(87) 



In this work, we used the replica cr-model technique to 
describe mesoscopic fluctuations and weak localization 
correction to the equilibrium supercurrent in Josephson 
junctions formed of a metallic wire between supercon- 
ducting leads. We considered two types of junctions: a 
chaotic dot coupled to superconductors by tunnel barriers 
(S-QD-S) and a diffusive wire (SNS) with transparent NS 
interfaces. In both cases we calculated the amplitude of 
supercurrent fluctuations and the weak localization cor- 
rection to the average current /(x) in different temper- 
ature regions at arbitrary ratios between A and -Bdwcii 
(given by Eg for an S-QD-S junction, and by Et for an 
SNS junction). 



For a quasi-one-dimensional SNS junction, we have 
found that mesoscopic corrections to the quasiclassi- 
cal Josephson current are "nearly universal" : SIc/Ic ^ 
—A^ffi^Ic/ Ic ^ Gq/G, where the exact coefficients in 
these relations are of the order of 1, being slightly de- 
pendent on the parameters of the junction. For a double- 
barrier S-QD-S junction, the weak localization correction 
vanishes, while mesoscopic fluctuations arc "less univer- 
sal": SIc/Ic ~ Gq/G for junctions with kTc ^ Eg, and 
is additionally suppressed for junctions with Eg ^ kTc- 

We also demonstrate that the approach of Ref. 5 sys- 
tematically underestimates the magnitude of mesoscopic 
fluctuations of the critical current by factors around 
2.5-2.8, both in the quasi-one-dimensional and three- 
dimensional cases. 

Theoretical predictions for the mesoscopic fluctuations 
of the critical current should be compared with experi- 
mental results. ^^'^"^ In Ref. 23, mesoscopic fluctuations 
of Ic were studied in a geometry of wide {W ^ L) long 
(Et ^ A) bar of a two-dimensional electron gas. Exper- 
imentally observed fluctuation magnitude (in the weakly 
localized regime I) is three times smaller than our re- 
sult (85). Experiment of Ref. 24 refers to the quasi- 
one-dimensional short (A ^ Et) wires, and should be 
compared with the result^^'^^ Sic — 0.30eA/ft. Again, 
theoretical prediction (^ 7 nA) appears to be several 
times larger than the experimentally measured magni- 
tude. This systematic discrepancy might be attributed 
either to the effect of electron-electron interaction or to 
non-ideal transparencies of the NS interfaces. 

Among the motivations to study mesoscopic fluctua- 
tions, it was suggested that they may induce the sign 
change of the critical current (vr-junction behavior).*^ 
However, such possibility was discarded in Ref. 43 if 
time-reversal symmetry is preserved in the normal part, 
and if interactions are absent. If the normal part is fer- 
romagnetic, time-reversal symmetry is broken. Then, 
with strong barriers between the wire and leads, it 
was shown that mesoscopic fluctuations are dominant 
in the Josephson relation compared to the quasiclassi- 
cal contribution.'*^ Within our approach this study could 
be easily reconsidered and extended by adding an ex- 
change energy term in the Hamiltonian of the wire. A 
strong Coulomb blockade in the central dot was also 
shown to provide a mechanism of 7r-junction formation 
via fluctuations.''^ We note that the replica cr-model is 
also appropriate for taking electron-electron interactions 
into account. This could be the subject of future study. 
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APPENDIX A: SOLUTION OF THE USADEL 
EQUATIONS FOR A WIRE 

1. Determination of the spectral current J and 6{0) 

For a fixed e > 0, dg and x, the values of the spectral 
current J and the Usadel angle Oq = 9{0) in the middle 
of the wire are determined from two equations: 



which can be solved as 



■ = 



r-^ dd 

le„ no)' 

1 d9 



(Ala) 
(Alb) 



where n{0) = [cos 6*0 - cos 6* + a(sin"^ 6^0 - sin"^ 9)]^/^ 
and a = ,P /Ae. These integrals can be converted to the 
standard elliptic integrals: 

_ 2F{^, k) 



\/a — c 



(A2a) 



(1 — a?')^a — c \/a — c 

' n(y>,fc^i^,fc) n(y,fc^i±|,fc) 

(l-a)(l-6) (l + a)(l + 6) 



(A2b) 



where 




Lp = arcsm 



{h — c){a — cos 9s) 



Note that our definition of elliptic integrals F{(p, k) and 
n((y9, n, k) coincides with that of Ref. 34. The same func- 
tions are often defined in a different way"'^: in Mathemat- 
ica, e.g., F{ip,k) = EllipticF[iy9, fc^]. and Il{(p,n,k) = 
EllipticPi[n, ip, fc^]. 

Equations (A2) should be solved numerically to obtain 
a and b (and hence J and 9q) for given e > 0, 9^ and 
X- The Usadel angles for e < can be obtained from 
Eqs. (49). 



2. Determination of 6{x) and ri{x) 

For each s we first have to find 9o and a as described 
above. Then for a spacial point Xi, the value 9i ~ 9{xi) 
can be found from the equation 



{a — c){b — cos 9i) 
{b — c){a — cos9i) 



(A4) 



cos 9j 



b~aY 



Y 



— ^ sn^{\xi\y^e{a - c), k), 
(A5) 

where sn(u, k) is the Jacobi elliptic function'^^. In Math- 
ematica, e.g., sn(M, fc) = JacobiSN[u, fc^]. 

Then we determine r]i ~ 'r]{xi) from Eq. (65): 



la 



d9 



So (cos 6*0 



no)' 



(A6) 



Converting this to elliptic integrals and using Eq. (A4) 
we get 



X, X, ^{b-a) Ili^^,k^^,k) 



(p-a) \xt\ y/a- c {p-a){p-b) 



(A7) 



where p = cos^o + a/sin'^6'o. Note that 7j{x) is an odd 
function of x. 

The functional determinants of A{x) involved in 
Eqs. (33) and (38) can be calculated numerically by intro- 
ducing a proper grid Xi. Then we discretize the Laplace 
operator in the operator A [Eq. (62)], and find amn, Pmn, 
rjm and rjn for each Xi, thus defining a finite matrix Aij. 
The functional determinants in Eqs. (33) and (38) can 
then be approximated by determinants of the matrix Aij 
which should be evaluated numerically. 



APPENDIX B: DERIVATION OF THE 
OPERATOR A FOR A WIRE 

Substituting Eq. (21) into (46) and expanding the ac- 
tion to the second order in W, we get: 



^(2 



G 



N 



1/2 



dxtr -iVW)^ + {J,Wr 

lOt^Q J-1/2 L 

- 2Jti[W,VW] - 2eUT3U''TiW^ 



(Bl) 



where J' = C/Vf/Vi. Using the decomposition (25) of 
W in terms of the d- and c-modes, we rewrite 5(2) in the 
form (61), where 

A„n{x) = -V^ + a„rn + A^„„/4 

iS2[A,,inV + (VAmn)/2] + (3m7i'^3 + 'Jmn'^l- (B2) 

Here, m = (e, a) and n = (e', b) are energy and replica 
indices, Sj=i,2,3 are the Pauli matrices in the {d, c)-space, 
amn is given by Eq. (63), and 



A 



(cos6'„V0„ 4- cos 6',„V (/)„), 



(B3a) 

= (sin0„V0„ sin6'™V0™ - V6l„V6'™)/2, (B3b) 
7mn = -(sin6'„V0„V6'm + sin6'mV<^mV6'„)/2. (B3c) 

Here, {9,n,4>m) and {9n,4>n) are the solutions of the Us- 
adel equations (48) at energies e and e', and phase dif- 
ference xi or X2, depending on the replica indices a and 
&, respectively. 
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In order to remove the first order derivative in 
Eq. (B2), we make a local unitary transformation in the 
(d,c)-space: (rf,c)^„j ^ Vjnn{d.,c)'^„, where 



■, , Cn "I" Cm --A ■ Cn ~^ C 
Vmn = COS t2^2 Sm 



Cm (a;) 



2 2 
ds cos 6'm(s)V0m(s). 



(B4a) 
(B4b) 



Such a rotation leaves detA and A^\^ invariant, while 
the operator A is transformed to 

A = 'V}rmAmnymm (B5) 

that can be written in the form (62) (tilde omitted) with 

V0„, 



r]m = — arctan — 



Cr, 



sin ^mV^, 

Taking the derivative of (B6), we come to Eq. (65) 



(B6) 



APPENDIX C: EQUIVALENCE OF EQ. (66) TO 
THE RESULTS OF REFS. 16 AND 17 



The result for the variance of the Josephson current, 
var/(x) = i5/^(x), obtained in Refs. 16 and 17 can be 
written in the form 



var/(x) 



1 /eA 



2 V Tih 



dkk{l 



-Tzk\ 



a{k) 



dx cos kx sin x 



cosliXA/cosh^ 2; — sin^ ^ 



|a(fc)P, (CI) 



(C2) 



In this Appendix we show that our expression (66) for 
((/(x)/(x'))) x' ~ X a-nd zero temperature gives the 
same result. 

We start with rewriting Eq. (66) in the form of integrals 
over C and C: 



x/2 



dC sinx 



dC sinx 



92 



sin C A / sin^ | - sin^ C sin C" W sin^ | - 



sin^ C 



; dCdC 



In 



sin(C - C") sin(C + C") 



C-C 



C + C 



Using Eqs. (2.29) and (3.6) from Ref. 16, we rewrite the second derivative of the logarithm as 



92 



■In 



sin(C - C) sin(C + C) 



-4 / dk 
/o 



k sinh kC sinh kC 



dCdC C-C C + C ./n e 

which allows to present Eq. (C3) in the form (CI) with a{k) instead of a{k): 

1 /"^/^ dC sinhfcC sinx 



h{k) 



sinh(7rfc/2) Jq 



sinCA/sin^ | - sin^ C 



(C3) 
(C4) 

(C5) 



The equivalence between our result and the result of Refs. 16 and 17 follows from the equality a{k) ~ a{k). To prove 
it, we extend integration in Eq. (C2) to the real axis, substituting cosfcx by e**""'^. Then we deform the integration 
contour to the upper half-plane and enclose all branch cuts of the square root: [i[Tm + tt/2 — x/2), i(7rn-|- 7r/2 + x/2)], 
n = 0, 1, . . . As a result, summation of £-("+1/2)7^/0 ^^^^ ^ yields sinh(7rfc/2) in the denominator, while integration 
along the branch cut reproduces the integral in Eq. (C5). Thus, a{k) = d{k) which establishes the equivalence between 
the two results. 



APPENDIX D: SOLUTION OF THE USADEL EQUATION FOR A WIRE CLOSE TO Tc 



Close to Tc, the solution of the Usadel equation (48) is given by Eq. (57). It corresponds to the following dependence 
of 6{x) and 4){x): 



smf 



sm W .T 



sinhtt; 



•\/cosh2Ka;(coshK — cosx) — (1 — cosh k cos x), 



tan 4>{x) 



tan ^ 

tanh KX 



tanh f 



(this form is valid both for e > and e < 0). The fimction 7]{x) defined in Eq. (65) is given by 

. tanh # 
tan ■ 



tan 77(2;) = — sign(e)-;— tanhrex. 



(Dl) 
(D2) 

(D3) 
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Rigidity of Gaussian fluctuations near the saddle point is determined by the operator A^^^.^ = — V^ + |£i| + |e2|+V^i62 j 
where the operator V can be expanded in the PauU matrices (73) with the coefficients 

V'^(x) — — (7—7 . — I (coshfti — cos Yi) cosh2Kia: — ^ I 7^ . — ) (cosh K2 — cos y2) cosh 2/t2a;, (D4) 
2 \|ei|suihKi/ 2 \|e2|smhK2/ 

V^(x) = —2- — r — — ■; — r — — f siKn(ei) cos — sin — sinh — cosh — sinh kix cosh K2X 
^ ' ei sinh €2 sinh K2 V ^^'^ 2 2 2 2 



+ sign(e2) sin ^ cos ^ cosh sinh cosh kix sinh K2X^ , (D5) 



V^(x) = 2-. — r — — ■; — r — — fsiu — siu — cosh — cosh — cosh kix cosh nox 
^ ^ |ei sinh Ki leal sinh K2 V 2 2 2 2 



sign(eie2) cos ^ cos ^ sinh sinh sinh kix sinh K2X ) . (D6) 



The supercurrent correlation function (71) involves 

= Tv{GVf = 2 dxdyGl,^{x,y) V\x)V\y). (D7) 

•^^1/2 i=0,l,3 

Since current fluctuations are determined by Reie2i it is convenient to symmetrize i? by introducing 

such that Y.ei.e2 -^^1^2 = Eei.ea ^eie2- 

Substituting Eqs. (D4)-(D6) into Eq. (D7), and taking the derivatives with respect to xi £^nd X2 we get 

92 /A Ki \ V A K2 ^ ^ 



^eie2 = 1 — f-^i 1 — r^T — T5iJ2sinxisinx2, (D9) 



dxidx2 \\f^i\smhKiJ \|e2|sinhK2^ 

resulting in Eq. (78). Here the function Te^^g is defined by the double integral: 

nl/2 



f r 

"^£1^2= / dx dy Gi:iE2i^Ty)^ cosh 2kix cosh 2^2^ 
J-1/2 L 



-1/2 

+ 2 cosh'^ cosh^ cosh cosh K2X cosh kij/ cosh Kay + 2 sinh^ sinh^ sinh kix sinh Ka^'^^ sinh Kiy sinh K2?/ 

— 2 cosh^ — !- sinh^ cosh kix sinh K2X cosh Kiy sinh K2y — 2 sinh^ cosh^ sinh kix cosh K2X sinh Kiy cosh K2y 

(DIO) 

In principle, these integrals can be calculated in a closed form. However, the resulting expression is too complicated 
and wc leave the integrals unevaluated. 
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